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ABSTRACT 



The morphologies of planetary nebulae are believed to be self-organized config- 
urations. These configurations are modeled by three-dimensional temporally self- 
similar magnetohydrodynamic solutions with radial flow, under the gravitational 
field of a central star of mass M. These solutions reproduce basic features, such as 
jets, point-symmetric knots, and filaments, through plasma pressure, mass density, 
and magnetic field lines. The time evolution function of the radial velocity starts 
as a slow wind and terminates as a fast wind. 



Subject headings: 
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1. Introduction 

The spherical observational features of planetary nebulae can be accounted for by the 
celebrated interacting wind model (Dyson and de Vries 1972, Kwok et al 1978), where the early 
slow wind of an intermediate mass star on the asymptotic giant branch phase is catched up 
by the following tenuous fast wind, generating the spherical features. The elliptic and bipolar 
features could also be reproduced by postulating the presence of a dense equatorial cloud (Kahn 
and West 1985, Mellema et al 1991), which has been imaged by high resolution instruments. 
Nevertheless, besides these spherical, elliptic, and bipolar features, images have also revealed 
point-symmetric knots (Miranda and Solf 1992, Lopez et al 1993, Balick et al 1993). These 
point-symmetric features plus the detection of magnetic fields in central stars (Jordan et al 
2005) call for a magnetohydrodynamic (MHD) approach of planetary nebulae (Pascoli 1993, 
Chevalier and Luo 1994, Garcia-Segura 1997, Bogovalov and Tsinganos 1999, Matt et al 2000, 
Gardiner and Frank 2001). Recently, Tsui (2008) has regarded the morphologies of planetary 
nebulae as self-organized structures, described by temporally self-similar MHD solutions in 
spherical coordinates with radial velocity flow for explosions. The connection between the 
physical process of self-organization (Hasegawa 1985, Zhu et al 1995, Yoshida and Mahajan 
2002, Kondoh et al 2004) and the mathematical analysis of self-similarity, using the Lagrangian 
radial label, has been discussed in detail. With axisymmetry, an equatorial plasma torus 
and a bipolar planetary nebula can be reproduced. In essence, current hydrodynamic and 
magnetofluid models attribute dynamic structures to the central star, or stars, as a cause to 
break the spherical isotropy, to account for the nebula morphologies, as a consequence. The 
self-similar approach relies on the global conservation properties of MHD plasma to reach self- 
organized configurations regardless of the initial conditions. The spherical isotropy of explosion 
is broken during the course of expansion because of the magnetic field, which is non-spherical 
by nature. 

We follow this self-similar view that considers astrophysical phenomena with ejection ori- 
gins as self-organized objects, whose configurations are solved by temporally self-similar MHD 
solutions (Low 1982a,b, Low 1984a,b, Osherovich et al 1993, 1995, Tsui and Tavares 2005, Tsui 
2006, Tsui et al 2006). This includes active galactic nucleus (AGN) jets (Tsui and Serbeto 
2007), classically treated as a steady state acretion-ejection MHD transport phenomenon from 
the acretion disk to the polar axis (Blandford and Payne 1982), and axisymmetric planetary 
nebulae. Here, we develope three-dimensional MHD solutions to describe additional planetary 
nebula features like jets (ansae), knots, and filaments. Unlike AGN jets that are extensive ob- 
jects in galactic scales and are collimated by some physical mechanism, such as plasma pressure 
driven collimation (Tsui and Serbeto 2007), nebula jets are much smaller objects and are not 
collimated. We remark that three-dimensional self-similar MHD solutions are very rare. To 
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our knowledge, there is only one published work which is devised specifically for interplanetary 
magnetic ropes (Gibson and Low 1998), which is not suitable for astrophysical phenomena. 
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2. Self-Similar Formulation 

The basic MHD equations in Eulerian fluid description are given by 



^ + V-(H = 0, (1) 

p{^ + (v-V)v} = JxB-Vp-p^f , (2) 

- = -Vx£ = Vx(i?xB), (3) 

Vx5 = /iJ, (4) 

V • £ = , (5) 

!-(-) + (v-V)(^-) = . (6) 



— * — * 

Here, p is the mass density, v is the bulk velocity, J is the current density, B is the magnetic 
field, p is the plasma pressure, \i is the free space permeability, 7 is the polytropic index, and 
M is the central mass which provides the gravitational field. 

We consider a radially expanding plasma and seek self-similar solutions in time where the 
time evolution is described by the dimensionless evolution function y(t). For this purpose, it is 
most convenient to think of Lagrangian fluid description, and consider the position vector of a 
given laminar flow fluid element r(t). Under self-similarity, the radial profile is time invariant 
in terms of the radial label 77 = r(t)/y(t), which has the dimension of r. Physically, r\ is the 
Lagrangian radial position of a fixed fluid element. With a finite plasma, the domain of rj is 
bounded by mass conservation 

< T]int < V < Vext ■ (7) 

As for the plasma velocity, we consider self-similar structures deriving from a spherically sym- 
metric radial velocity which can be written as 



(8) 
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Our self-similar parameter 77, defined through the Lagrangian fluid label, explicitly represents 
the fluid velocity by the time evolution function y(t). This evolution function will be solved 
self-consistently with respect to the spatial structures of the plasma. We emphasize that self- 
similarity, as a method, can be applied in different ways other than the one we use here. For 
example, Lou and his collaborators have treated an aggregating fluid under its self-gravitational 
field with a similarity variable x = r(t)/at, where a is the sound speed, for isothermal fluid, 
7 = 1, (Lou and Shen 2004, Bian and Lou 2005) and for a polytropic gas, 7 > 1, (Lou and 
Wang 2006, Lou and Gao 2006) to study relevant astrophysical phenomena. Extensions to 
a magnetofluid have been considered by Yu and Lou (2005) and by Lou and Wang (2007). 
Because of the linear dependence on time, this similarity variable x refers to a reference frame 
moving at speed a, which is different from the radial plasma flow velocity v. Furthermore, 
different from our similarity variable 77, x here is not the Lagrangian label of a given fluid 
element. For this reason, the convective derivative remains explicit in the x representation. As 
a result, the similarity variable of Lou amounts to finding the plasma structures in an adequate 
moving frame in the Eulerian x fluid description. This resembles the analytic technique of 
going to a moving frame to look for stationary profile solutions for nonlinear phenomena such 
as nonlinear Alfven waves, solitons, etc. Because of this fundamentally different definition of 
self-similarity, the nature of the phenomena intended to describe is different. Our Lagrangian 
label self-similarity parameter is aimed to find spatial plasma configurations, and to determine 
the radial plasma flow velocity, consistent to the spatial configurations, through the evolution 
function. Since we are considering an isotropic radial plasma flow, a natural solution would 
be a hydrodynamic one-dimensional expanding plasma, with radially dependent mass density 

— * — * 

p and plasma pressure p, and with B = and J = 0. Nevertheless, this one-dimensional 
hydrodynamic solution is highly unlikely because magnetic field fluctuations can be generated 
from current density fluctuations, even in the absence of a pre-existing magnetic field. With 
the magnetic fields, which are basically a two- or three-dimensional structure, coupling to the 
plasma will generate likewise two- or three-dimensional p and p. 

The independent variables are now transformed from (r,6,<f),t) to (i],6,(f),y). We now 
proceed to determine the explicit dependence of y in each one of the physical variables with 
this radial velocity using functional analysis. First, making use of Eq.(8), Eq.(l) renders 

! + is (r v) = (t+A + 4^) 

at r 2 or at Or or r 



(9a) 
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To reach the second equality, we note that the first bracket in the first equality corresponds 
to the total time derivative of an Eulerian fluid element which amounts to the time deriva- 
tive of a Lagrangian fluid element. As for the second bracket, it can be reduced by using 
v = dr/dt = rjdy/dt and dv/dr = (l/y)(dy/dt). Solving this equation for the y scaling by 
separating the time part gives 

p(f,t) = -^pM,0) • (96) 

y 

As for Eq.(6), with a$F = Qo/p 7 ) where «o is a constant that carries the physical dimension so 
that F is a dimensionless function, it follows 

OF OF dFdy , . 

(4) = <* F(r,t) = -^a FM>) ■ (106) 

As for Eq.(3), with the aid of Eq.(5), the magnetic fields are 

dB r 1 d . ,„ . ,dB r dB rx 2v „ 

= m. dy + 25,^ = SB, + 2B, ^ = 
ay at y dt dy y dt 

B r (f,t) = \B r (ri,6,<l>) , (H6) 

y 

dt r dr dt dr r dr 



dB e dy + 25^ = + = (12a) 

dy dt y dt dy y dt 



B e (f,t) = ^B e (v,0A) , (126) 



^ + ~(^) = 0, (13a) 
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B^t) = -^(77, 0,0) . (136) 

y 

Making use of Eq.(9b), we derive the plasma pressure from Eq.(lOb) 

p(f,t) = 4^M, M) = = -^pM, ■ (14) 

y y y 

Since Eq.(12a) and Eq.(13a) are of the same form, we conclude that, under self-similarity, B^ 
is a linear functional of Bg with 

= kB e . (15) 

Making use of Eq.(4) to eliminate the current density in Eq.(2), we get the momentum equation 
which has three components. The 0, 9, and r components are respectively 

S4-(^sin^)-^]-S r [^-sin^-(^)]- y 4 -^-(^) = , (16) 



fc^(^sin0)-^]+3 r [sm^ = , (17) 



, - r 1 dB r 19.,-.. -.19. - , 1 dB r , , ~ 9 , , 
hB e \——^ - - — (r,kB e )] - B e [- — ( V B e ) - -— r -] - y^—^p) 
r)sm9 90 r)Or) r) Or) r) oO Or) 

_ 2 d 2 y _GM 

= wy^i + w-r- ( 18 ) 

We have reduced the general set of time- dependent ideal MHD equations, Eqs.(l-6), to a set of 
self-similar equations with appropriate time scalings, Eqs.(7-13). The general ideal MHD set 
has nonlinear terms of convective type (v ■ V). By using the fluid label description, the (v ■ V) 
convective terms are absorbed in the Lagrangian time derivative representation. The structure 
of the nonlinear terms, absorbed in the Lagrangian fluid label formulation, will appear in the 
r) profile of the system. 
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3. Jet Structures 

After this self-similar formulation, we have to solve Eqs. (16-18) for the self-similar con- 
figurations. To solve these equations, we first separate the radial variable from the other two 
variables by writing 

B r {ri,e,(j>) = AoR{rj)B r {0,<l>) , (19) 

and likewise for Bg, where B^ = kBg. Furthermore, we take 

p(ri,e,<f>) = PoR 2 (v)P(0A) ■ (20) 

Here, R(rj), B r (9,4>), and p(6,<f)) are dimensionless functions, and A and po carry the dimen- 
sions of magnetic field and pressure respectively. We have taken an R 2 (r]) dependence because 
plasma pressure is a quadratic positive quantity. We take A — 1 for unit amplitude magnetic 
fields, such that po is relative to this amplitude. Specifically, we take 

R(rj) = (a V y n , (21) 

to represent a power law decaying field with distance, where a is a normalizing parameterof rj. 
Considering (4 — 37) = 0, Eqs. (16,17) are respectively 

d ~ OB ~ 1 

B e [— (kB 9 smO)--^]-B r (n-l)sm6kB g -— (-B 2 r + fip p) = , (22) 



kB e [— {kBo sin6) --^]+B r {n-l)sm6B e + smO— {-B 2 + fip p) = , (23) 



We now separate the azimuthal dependence by writing 

fl r (M) = ®r(0)H(f>) , (24a) 

B e (e,<f>) = Q(0)*(<j>) , (246) 



p(M) = ej(0)$ 2 (0 , 



(24c) 
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where = kBg. These functional dependences give Eq.(22) as 

8 1 <9$ 

e— (fee sine) _e r (n-i) sin eke = [e 2 - (e r 2 + 2 m e;)]-— 



Considering the separation constant im, such that 



^ ////'!» , (25a) 



$(0) = e +tm4) = cosm0 + 2sinm0 , (256) 



we then have 

e— (e sine)- e r (n-i) sin = +^m[e 2 - (e 2 + 2 m ©;)] , (26) 

which identifies fc = +i. Following the same procedures, Eq.(23) reads 

6 — (0 sin^ -6 r (n- 1) sin #6 = +m6 2 + - sine— (0 2 + 2/xp Qp) • (27) 
cw 2 ae 

The left sides of these two equations are the same, which allows the right sides be equated to give 

<9G 2 

sine^l = _ 2m 9 2 , (28a) 



e 2 = (e 2 + 2fx Po e 2 p ) . (286) 

To solve for 6 2 , we integrate Eq.(28a) to get 

ln0 2 = -2m / — — = -2m / — Y - = +ln m , 29 

7 sine J sin 2 e v l-aT 

where we have multiplied and divided the right side by sin 6 to implement the integration, and 
x = cose. This solution of 6 2 is singular at x — +1 for positive m, and x — — 1 for negative 
m. Such solution gives jet features on magnetic field lines and plasma density. 
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4. Special n=2 Case 

To get 6, instead of solving either Eq.(26) or Eq.(27), we make use of Eq.(5) to get 

V-B = ^(i^)){-(n-2)4 + ^[^(sin^) + A(^)]} = , (30a) 

d 

— (sin#0) = mQ + (n-2)sin#0 r . (306) 

Subsituting Eq.(30b) into Eq.(26) and Eq.(27) respectively gives 

sin#ee r = mQl , (31a) 



200 r = -^0* . (316) 

It can be shown readily that these two equations are consistent to Eq.(28a). We note that 
Eq.(30b) couples with r . To decouple these two functions, we consider the special, but 
probably practical, case of 

n = 2 . (32) 



Multiplying over by sin 6 and defining P{6) = sin 6*0, Eq.(30b) can be integrated to give 

P(x) = ([^) m/2 , (33a) 
(I - x^™- 1 )/ 2 

«<*> = jTTTfi^ • < 336) 

e,.(x) = m(ii^) 3m/2 , (33c) 

X Ob 

2WoO;(x) = [(^) m - m\\±^f m ] > , (33d) 

X Ob X Ob 

where Eq.(33c) is obtained from Eq.(31a), and 9p(x) can be recovered from 0*(x) of Eq.(28b). 
From these solutions, we see that B r and p are singular at x — +1 through Q r (x) and ®p(x), 
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whereas Bg and are singular at x = —1 through Q(x). We require the singularities be 
integrable in x, which demands the power of the singularities be less than unity. By inspection 
of the terms, we conclude that m < 1/3. Let us take 

m = \ , (34) 

to show Ql(x) in Fig.l in a polar plot, Q r (x) and Q(x) in Fig. 2 and Fig. 3 respectively. Figures 
1 and 2 show an integrable singularity at x — +1, or 9 = 0. Furthermore, Fig.3 shows a 
weaker singularity at x — +1 than the one of Fig. 2. However, Fig.3 shows another singularity 
at x — — 1, while Fig. 2 is regular at that location. As for 9p(x), the second term of Eq.(33d) 
exceeds the first term when 

m(l±^r > 1 , 

1 — x 

or as x gets very close to unity. Writing x — (1 — 5), we get (5/2) m < m, or (5/2) < m y l m . 
With Eq.(34), we have 5 = (1/2) 7 . Consequently, Q 2 {x) gets smaller as x approaches unity, 
and it vanishes at x — (1 — 5), beyond this point, it gets negative. To interpret this negative 
plasma pressure amplitude, Qp(x) now points backwards to the x — — 1 direction in a polar 
plot, giving a jet structure as in Fig.4. This jet structure is along the radial magnetic field 
given by Q r (x). 

As for the radial component, Eq.(18), with (4 — 37) = 0, it reads 

2 d 2 y GM .l^oxfsr, 1 9B r dB ri n A _ 
at 1 i] 2 i] smO o<p oO r\ 

= {-R 2 )B e \-mQ r + sin^]$ + 2n(-i? 2 ) m (0 p $) 2 . (35) 

With a as the separation constant, we have 

(-R 2 ){B e [-mQ r + sin0^]$ + 2n m (e p $) 2 } = npi^- + a V ) , (36a) 

7] OU T] z 



d 2 y a 



(366) 



We write the mass density as 
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P(V,0,<P) = p K(vWA) , (37) 

where po carries the dimension of mass density, and R(^) and p(8, 0) are dimensionless func- 
tions. We can identify immediately from Eq.(35a) that 

R 2 

R(V) = 777TT7 n\ , (38a) 

v " (GM/rj + at] 2 ) v ' 



ppop = {B e [-m& r + sin^^]$ + 2n/ip o (0 P $) 2 } = [2npp & 2 p - 2m9e r ]$ 2 . (386) 

From the first of these two equations, we deduce that a is positive, such that R(^) is analytic. 
As for the second equation, it reads 

1 _|_ r 1 _|_ r ( 1 _|_ r )(2m-l)/2 

= {»Kf^r - ™ 2 (^) 3 "'] - .%o^ m ^ > ■ < 39 > 

This mass density has a singularity at x = +1. With m = 1/4, we note that the numerator 
of the last term in Eq.(38), (1 + a; )( 2m " 1 )/ 2 ; has a negative power. This gives a singularity at 
x = —1. Since the power of this singularity is —1/4, it is also integrable as well. The mass 
density distribution is shown in Fig. 5 with also a jet structure in the x — — 1 direction, because 
of the radial magnetic field. If we consider m < negative, we should point out that Figs. (1-5) 
would turn upside down. Consequently, the jet structures would be on both sides of the polar 
axis. Some examples of these jets appear in M2-9 Twin Jet Nebula, CRL 2688 Egg Nebula, 
NGC 3242, NGC 6826, NGC 7009. As for the time evolution function of Eq.(35b), we multiply 
over by dy/dt to get the first integral as 

(§) 2 = 2(tf~)>0, (40a) 

where H is an integration constant. Knowing that a is positive, and y(0) = 1 by definition of 
Lagrangian fluid label, and the right side of this equation has to be positive, we conclude that 
H > a > 0. This gives the plasma wind velocity a slow start initially with (dy/dt) 2 = 2(H — a) 
which evolves into a fast terminal wind of 



(406) 
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5. Knot and Filament Structures 

With the momentum equation solved, the magnetic fields are given by 

B r = +R(r])Q r cosm0 , (41a) 
B = +R(ri)Qcosm<j) , (416) 
B$ = -R(rj)Q sin m0 . (41c) 



The magnetic field lines are given by 

G r cos m0 cos 7770 G sin 777<; 



dr] r]d6 77 sin 6d(f) ' 

which can be written as 



(42a) 



dn (l+x) 2m , l + :r 2m+1 d0 , in ,. 
-L = -m) L—dx = -m) '- — — . 426 

7/ (1-X) 2m (1 - X) 2m - 1 tan 7770 



With 7/0 as an integration constant, the first equality gives 

In ^ = -m ! [ l + X yT dx , (43a) 

7/o J (1-X) 2m 



which can be integrated numerically, as is shown in Fig. 6. Since the singularity is integrable, r\ 
is finite at x — +1. The second equality corresponds to 

d6 d0 



sin 9 tan m0 ' 

and can be integrated to give 

|7±-^ = ^(smm0) 2 , (436) 

where K is an integration constant. With the interval of x between (— 1,+1), or 9 between 
(it, 0), the left side of Eq.(43b), (1 + x) m /(l — x) m labelled on the left axis, covers an interval 
(0, oo) with positive m, and is plotted in Fig. 7 against x between (— 1,+1) labelled on the 
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bottom axis. The right side, K(smm<f)) 2 labelled on the right axis, is also plotted against m0 
on the top axis, with a scale between and tt/2, in the same figure. In order to view the 
mapping between the right side and the left side, we assign a large constant K. With m — 1/4 
and K = 3, as x departs from — 1, or from ir, mcj) departs from 0, and it maps a root of <p. As 
x approaches +1, or 9 approaches 0, m0 reaches tt/2, which takes over (0, 2tt). On the return 
path of the field lines, x decreases from +1 back to —1, bringing mcj) from 7r/2 to n along the 
descending branch of K(smm<f)) 2 , which takes over 2tt, 4tt. This descending branch, which is 
the continuation of Fig. 7, is not shown. To summarize, the field lines starting at x — — 1 and 
i] = +1 go through the x = (—1, +1, —1) cycle once, with i] = (+1, +2.3, +1), while completing 
the m4> = (0, 7r) cycle once, covering <fi = (0, 4ir), before closing on themselves again. This 
generates helical field lines, as shown in Fig.8, on the surface of revolution of Fig. 6. Since the 
singularities at x — ±1 are integrable, and also because of the circulating nature of the fields 
Bq and 13$, the magnetic field lines converge to 77 = 1 at x — — 1 axis and to 77 = 2.3 at x — +1 
shown in Fig. 6. These locations correspond to point-symmetric magnetic knots, where 
the field strength is infinite. As for the filaments, they correspond to the helical magnetic field 
lines in space, as shown in Fig.8. These same field lines give a different shape when they are 
viewed at different orientations, such as in Fig. 9 and Fig. 10. Further field lines can be generated 
with, for example, m = 1/3.5 = 2/7 = 4/14. In this case, the field lines starting at x — — 1 and 
77 = +1 go through the x = (—1, +1, —1) cycle once, with rj = (+1, +2.3, +1), while completing 
the m4> = (0, 7r/2, it) cycle once, covering <fi = (0, 7tt/4, 7n/2). Since the lowest 2-7T multiple of 
the <p cycle is four, the field lines have to complete four cycles of x and of m(p to make covering 
(0, 7n, 147r), such that the field lines can close on themselves again. With different values of 770 
of Eq.(43a), we can fill up the space with shells of field lines of some m less than 1/3 up to r) ext 
of Eq.(7), and dot the polar axis with a line of magnetic knots. Despite the radial component 
at x — +1 giving the plasma jet structures, the field lines are closed at rj = 2.3 because of 
the circulating nature of the meridian and azimuthal components. With m < negative, the 
mirror images of the knots and field lines can be superimposed on those with m > positive, 
generating lines of knots and concentric shells of magnetic field lines. Examples of knots can 
be found in M2-9 Twin Jet Nebula, NGC 5307, and filaments in MyCn 18 Hourglass Nebula, 
NGC 6543 Cat's Eye Nebula, NGC 2392 Eskimo Nebula, M2-9 Twin Jet Nebula, NGC 6543. 
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6. General n Case 

We now solve Eq.(30b) for an arbitrary n. Combining Eq.(30b) an Eq.(31a), we get 

(1 -x 2 )^-P 2 {x) + 2mP 2 (x) = -2m(n - 2)(1 - x 2 )Ql(x) , (44a) 

where P(x) = (1 — x 2 ) l l 2 Q(x) and Q 2 (x) is given by Eq.(29). If the boundary condition is 
known, this equation can be integrated numerically to get 



P 2 (x) = -2mJj-0^ + (n-2)Ql(x)}dx + P 2 (-l) . 



(446) 



To obtain the boundary condition at x — — 1, we note that the right side of Eq.(44a) is equal 
to -2m(n — 2)(1 + x) m+1 /(l — rr)" 1 " 1 . Since m < 1 positive, this term vanishes at x — — 1 
and at x — +1 as well. As a result, in the neighborhood of x — ±1, P 2 (x) is described by the 
homogeneous version of Eq.(44a), where the right side is null. The homogeneous solution can 
be solved readily as 

P\x) = (|^r , (45) 

which provides the needed boundary condition. The self-similar functions are, therefore, given 
by 

9(x) = , (46a) 

= ^ , (466) 

2 m e 2 p (x) = [el(x) - e 2 r (x)\ > o , (46c) 



As for the magnetic field lines, they are described by 



In 1. = _ m f |^4da; , (47a) 



(476) 



Vo J P 2 (x) 

de _ d<j) 

sin 6 tan m(f> 
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The second equation remains the same, while the first equation can be integrated numerically 
to get the field lines. With n — 3, the corresponding self-similar functions of Figs. 2-6 are 
evaluated anew, and are presented in Figs. 11-15 respectively. From Fig. 12 and Fig. 15, we can 
see that the source term in Eq.(44a) makes Q(x) and the (rj — x) mapping more symmetric. 
The corresponding magnetic field lines are also shown in Figs. 16-18, with knots, x — ±1 and 
r\ — +1, at equal distance from the center. 

To conclude, we have regarded the morphologies of planetary nebulae as self-organized 
configurations. These configurations are modeled by temporally self-similar MHD solutions. 
To complement an earlier publication (Tsui 2008) for axisymmetric features, we have presented 
a three-dimensional self-similar model with 7 = 4/3, which reproduces features like jets, point- 
symmetric knots, and filaments, through plasma pressure, mass density, and magnetic field 
lines. The time evolution function of the self-similar solutions starts the plasma expansions 
as a slow wind, and terminates as a fast wind. With this three-dimensional model, which 
completes the earlier two-dimensional axisymmetric model, we have covered most of the existing 
features in planetary nebulae with this self-similar approach for self-organized configurations. 
Considering that the extragalactic AGN polar jets could be accounted for on the same ejection 
basis (Tsui and Serbeto 2007), through a plasma pressure driven collimation process, which 
differs from the classical accretion-ejection spatially self-similar steady state MHD transport 
model (Blandford and Payne 1982), we believe temporally self-similar MHD configurations are 
universal manifestations of self-organized astrophysical ejection phenomena. 
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Fig. 1. — The function Q*(x) with x = cos 9 is shown in a polar plotted indicating an integrable 
singularity at x = +1 for a jet structure. 
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Fig. 2. — The function Q r (x) with x = cos 9 is shown in a polar plotted indicating an integrable 
singularity at x = +1 for a jet structure. 
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Fig. 3. — The function O(x) with x = cos# is shown in a polar plotted indicating integrable 
singularities at x — +1 and at x = — 1 for knot structures. 




Fig. 4. — The function Qp(x) with x = cos 9 is shown in a polar plotted indicating a plasma 
pressure jet in a very narrow cone about x — +1. 
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Fig. 5. — The function p(x) with x = cos 8 is shown in a polar plotted indicating a mass density 
jet in a very narrow cone about x — +1. 
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Fig. 6. — The magnetic field line r] — 6 dependence is shown in a polar plotted. 




Fig. 7.- 



The magnetic field line — 6 dependence is shown in a parametric plotted. 
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Fig. 8. — The three-dimensional magnetic field lines, wounded on the surface of revolution of 
Fig. 6, are viewed parallel to the x-y plane at about 45 degrees. 




Fig. 9. — The three-dimensional magnetic field lines, wounded on the surface of revolution of 
Fig. 6, are viewed down the z axis. 
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Fig. 10. — The three-dimensional magnetic field lines, wounded on the surface of revolution of 
Fig. 6, are viewed along the x axis. 
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Fig. 11. — The function Q r (x), with n = 3, is shown in a polar plotted indicating an integrable 
singularity at x = +1 for a jet structure. 
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Fig. 12. — The function Q(x), with n = 3, is shown in a polar plotted indicating a more 
symmetric structure. 
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Fig. 13. — The function Qp(x), with n = 3, is shown in a polar plotted indicating a plasma 
pressure jet in a very narrow cone about x — +1. 
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Fig. 14. — The function p(x) with n = 3 is shown in a polar plotted indicating a mass density 
jet in a very narrow cone about x — +1. 
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Fig. 15. — The magnetic field line rj — 6 dependence, with n = 3, is shown in a polar plotted 
showing a symmetric structure. 
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Fig. 16. — The three-dimensional magnetic field lines, with n — 3, are viewed parallel to the 
x-y plane at about 45 degrees. 




Fig. 17. — The three-dimensional magnetic field lines, with n 



= 3, are viewed down the z axis. 
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Fig. 18. — The three-dimensional magnetic field lines, with n = 3, are viewed along the x axis. 



